Nuclear Stopping and Sideward Flow Correlation from 0.35v4 to 200^4 GeV 



OO 
O 

o 

(N 

^■ 
CD 

m 



> 

m 

G\ 
O 
00 

O 



Xiao-Feng Luo/'0 Ming Shao/ Xin Dong,^ and Cheng Li^ 

^ University of Science and Technology of China, Hefei , Anhui, 230026, China 
(Dated: September 13, 2008) 

The correlation between the nuclear stopping and the scale invariant nucleon sideward flow from 
SIS/GIS to SPS/CERN energies is studied within Ultra-relativistic Quantum Molecular Dynamics 
(UrQMD). The universal behavior of the two experimental observables for various colliding systems 
and scale impact parameters are found to be highly correlated with each other. As there is no phase 
transition mechanism involved in the UrQMD, the correlation may be broken down by the sudden 
change of the bulk properties of the nuclear matter, such as the formation of Quark-gluon plasma 
(QGP), which can be employed as a QGP phase transition signal in high energy heavy ion collisions. 
Furthermore, we also point out that the appearance of breaking down of the correlation may be a 
powerful tool to search for the critical point on the QCD phase diagram. 

PACS numbers: 25.75.Ld, 25.70.Pq 



In the recent years, the main aim of ultra-relativistic 
high energy heavy ion collisions (HICs) performed at 
SPS/CERN {^/s^ - 10 A GeV) and RHIC/BNL 
(v^snn ~ 200A GeV) is to search a new form of matter 
with partonic degrees of freedom, the so-called Quark- 
gluon plasma (QGP) [E 0, i, Q. Although great ef- 
forts have been made, no dramatic changes of experi- 
mental observables, such as jet-quenching, elliptic flow 
and strangeness enhancement, have been observed yet 
and it is hard to make a solid conclusion for the hap- 
pening of the QGP phase transition Recently, the 
energy scan program is proposed for RHIC /BNL to per- 
form HICs experiments with lower cm. energy to search 
for the critical point [1, 0, Q > which is as a endpoint of 
the first order phase transition line on the QCD phase 
diagram. If the critical point exists, it should appear on 
the QGP transition boundary at higher baryon chemi- 
cal potential and lower colliding energy M- To ex- 
tract the QGP phase transition signal, a large amount 
of possible experimental probes, such as particle ratio 
and collective flow etc. have been proposed. The time 
evolution of temperature and baryon chemical potential 
of the different colliding nuclei with various colliding en- 
ergy would be mapping much broader T — hb region on 
the QCD phase diagram than a single nuclei. However, 
it is complicated to uniformly and systematically obtain 
unambiguous experimental signal for QGP phase transi- 
tion and also mark the location of the critical point on 
the QCD phase diagram, one of the possible choice is to 
keep insight in universal correlation pattern of two exper- 
imental observables for various colliding systems (system 
size and beam energy). Thus, the complication of vari- 
ous colliding systems dependence of searching the phase 
transition signals can be reduced. 

In this work, the correlation between the nuclear stop- 
ping and scale invariant nucleon sideward flow within 
the framework of UrQMD model from SIS/GSI to 



SPS/CERN energies for various scale impact param- 
eters < 6o = b/bmax < 1 have been found, not 
just for global fixed impact parameters as in Ref . [T]| . 
The scale invariant nucleon sideward flow is defined as 
F {bo) = d{< p-JA > /pPD/5(2/c™/y?r , first 

proposed in Ref. the linear fitting slope of the nor- 
malized rapidity dependence of the normalized average in 
reaction plane transverse momentum with fitting range 
of -1< ycTn/Vcm^ <lj where < Pcm/A > is the average 
transverse momentum projected in the reaction plane per 
nucleon and p^^-' is projectile momentum in the center 
of mass system (c.m.s.), j/cm and y^™-^ are the nucleon 
rapidity and projectile rapidity in c.m.s., respectively. 
The nuclear stopping ratio i? as a measurement of de- 
gree of stopping of colliding nuclei at scale impact pa- 
rameter bo, suggested in Ref. [H, [l^, is expressed as 
Ribo) = f \Pti\/J2i bzili where pu and p^i are trans- 
verse and longitudinal momentum of the ith outgoing 
particle in the c.m.s., respectively. A colliding system 
dependence variable /9,„f, is also introduced to be a nor- 
malization factor for later calculations. It is defined as 
p,nb{A,Eiab) = MB{Q)*ul';^^ /A''/^ where the MB(0) 
stands for the meson to baryon ratio for central collision 
(6o = 0), uP™^' = pP^hl^^ is the spatial component 
of four-velocity of the projectile in the c.m.s and A is 
the mass number of a nuclei in the symmetric colliding 
system. The overlapping volume of two collide nuclei 
and the nuclear passing time for central collisions respec- 
tively satisfy V (x A^ and t 
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where r is the radius of nuclei. Thus, we have pmb oc 
MB{0)/{V * tpass), standing for the meson to baryon 
ratio per unit volume, per passing time in the central 
collisions, which is used to characterize the strength of 
particle production at early stage 15]. 

The UrQMD model [lB| used here is a type of nu- 
merical transport models, which is based on the quark, 
di-quark, string and hadronic degrees of freedom. It in- 
cludes 50 different baryon species (nucleon, hyperon and 
their resonances up to 2.11 GeV) and 25 different me- 
son species. Two types of equation of state, the hard 
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FIG. 1: (Color Online) Left panels: The excitation function 
of the central nuclear stopping ratio, the variable pmb and 
the semi-central (0.3 < feo < 0.4) scale invariant nucleon side- 
ward flow. Right panels: System size dependence of the three 
experimental observables. 



EoS with incompressibility K = 380MeV(Only for beam 
energy up to AA GeV) and cascade are contained in the 
UrQMD model. The model has successfully been applied 
to reproduce the experimental results from SIS/GSI to 
SPS/CERN energies [13 . 

A group of symmetric colliding nuclei with five 
pairs: i^^Au+i^^Au, i29xe+i29Xe, 96Ru+96Ru, s^Ni+ssNi, 
'"'Ca-|-'*tla are combined with 30 and 12 incident kinetic 
energies, respectively. The first combination including 
total 150 = 5x 30 colliding systems with 30 beam energies 
per nucleon from 0.35 to 200 GeV (0.35, 0.5, 0.66, 0.83, 
1.0, 1.5, 2.0, 3.0, 4.0, 5.3, 6.6, 8.0, 10.0, 10.93, 11.9, 12.9, 
13.93, 15, 16.9, 17.92, 18.95, 20, 24.22, 36.0, 55.0, 76.92, 
102.33, 131.33, 163.36, 200.0) for each pair of the collid- 
ing nuclei and the second one including total 60 = 5 x 12 
colliding systems with 12 beam energies per nucleon from 
0.35 to 3.9 GeV (0.35, 0.5, 0.66, 0.83, 1.0, 1.5, 2.0, 2.35, 
2.7, 3.1, 3.5, 3.9) are researched with cascade and hard 
EoS of UrQMD, respectively . 

Fig. 1 shows the beam energy and system size depen- 
dence of central (&o = 0) nuclear stopping ratio i?(0) 
and the predefined variable pmb as well as the semi- 
central (0.3 < 6o < 0.4) scale invariant nucleon sideward 

flow F (&o) with the cascade EoS. In the left panels of 

Fig. 1, the i?(0) and F (bo) both decrease monotonously 
with the beam energy per nucleon from 0.35 to 200 GeV 
for three pairs of symmetric colliding nuclei: Au-|-Au, 
Ru-f Ru and Ca-|-Ca, and larger nuclear stopping ratio is 
observed for heavier colliding nuclei than the lighter one 
for a fixed beam energy. More detailed information on 

the system size dependence of -R(O) and F (bo) are illus- 



FIG. 2: (Color Online) The correlation between the central 
nuclear stopping ratio and scale invariant nucleon sideward 
flow, with bo varying from 0.2 to 0.6 and an interval of 0.1, 
are calculated for the mentioned 150 colliding systems within 
the cascade EoS. The solid line in each panel is the linear fit 
of the corresponding correlation. 

trated in right panels of Fig. 1. Both i?(0) and F (bo) 
increase monotonously with mass number A, which is 
proportional to the size of the colliding system. The 
defined variable Pmb, increasing with beam energy and 
decreasing with system size, is also shown in lower panel 
of the Fig. 1. 

With cascade EoS in UrQMD, the observables: cen- 
tral (6o = 0) nuclear stopping ratio i?(0), pmb and scale 

invariant nucleon sideward fiow F {bo), with non-zero bo 
satisfying < bo < 0.8 and with an interval of 0.1, are 
calculated for the mentioned 150 colliding systems with 
beam energy per nucleon from 0.35 to 200 GeV. After the 
logarithmic operations are performed on both normalized 
nuclear stopping R{Q)/pmb and normalized scale invari- 
ant nucleon sideward flow F (bo)/ Pmb, the resulted two 

variables, ln[R{0)/pmb] and ln[F (bo)/ Pmb] show strong 
universal correlation for various colliding systems with a 
given bo bin. For illustration, the correlation with the bo 
from 0.2 to 0.6 and the corresponding linear fit line are 
shown in Fig. 2. For the hard EoS case, the results of 
mentioned 60 colliding systems with beam energy per nu- 
cleon from 0.35 to 3.9 GeV are also shown in Fig. 3. The 
superposed solid line shown in Fig. 2 and Fig. 3 are the 
results of the linear fit for the corresponding correlation. 

By the linear fit of the correlation in Fig. 2 and Fig. 
3, the analytic relation between the two variables: the 

i?(0) and F {bo) can be expressed as: 

, r^(0)i . , J {bo), , 

ln[ — —\ = L X ln[ — ^ — -] + m (1) 

Pmb Pmb 
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Hard EOS: 0.35A to 3.9A GeV 
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FIG. 3: (Color Online) The correlation between the central 
nuclear stopping ratio and scale invariant nucleon sideward 
flow, with 6o varying from 0.2 to 0.6 and an interval of 0.1, 
are calculated for the mentioned 60 colliding systems within 
the hard EoS. The solid line in each panel is the linear fit of 
the corresponding correlation. 

, where the two fitting parameters L and m are intro- 
duced to represent the slope and intercept, respectively. 
Generally speaking, there is nothing particular for any of 

the two variables R and and they are of equal impor- 
tance. Actually, it is found that not only central stopping 
ratio, but also non-central nuclear stopping ratio is corre- 
lated with the scale invariant nucleon flow F (bo), which 
means for two independent scale impact parameters &^ 
and Bq , the relation between the corresponding 

and F (^o ) can be expressed as: 

ln[«]^Lxln[M)H.. (2) 

Pmb Pmb 

The fitting parameters L and m, depending on both 
and b^ for the cascade and hard nuclear EoS cases, 
are illustrated in Fig. 4. In the two upper panels of 
Fig. 4, the L shows almost no dependence on the b^ 
and bg and it is larger for hard EoS than the cascade 
one. Thus, it can be regarded as a constant parame- 
ter to characterize the nuclear EoS. The 6^ and 6^ de- 
pendence of m are also illustrated in the two lower pan- 
els of the Fig. 4, which is strongly affected by Bq , Bq 
and also by the different nuclear EoS. The two param- 
eters are both colliding systems independence, as they 
are both the universal fitting parameters for various col- 
liding systems. The parameter m = to(6^,6^) is only 
as a function of Bq and 6q', and the so-called correla- 
tion function C{Bq,Bq) = e~"^'^''o '''o^ is defined to de- 
scribe the correlation strength between the nuclear stop- 
ping ratio R{Bq ) and scale invariant nucleon sideward 



FIG. 4; (Color Online) The dependence of fitting parameters 
L (slope: two upper panels) and m (intercept: two lower 
panels) on scale impact parameters and for cascade 
and hard nuclear EoS. 

fiow F {b^)- Two new variables: i?*(6^) = R{b^)/pm,b 

and = F {1)^)/ p^^b are defined for the simplifi- 

cation of equ. (2). Then, it can be rewritten as: 

F''{b^)=(R*{b^)C{B^,X)Y 

, where the correlation function, < C{Bq,Bq) < 1, is 
only related to the Bq and Bq for given nuclear EoS. 

The colliding system as well as scale impact parameter 
dependence of single experimental observable are further 
investigated for any specific colliding system. From the 
equ. (3), for any given Bq and 6^, the two correlative 
observables are respectively calculated with two different 
scale impact parameters, {Bq^,Bq'^) and (5q^^,6^^), then 
we have: 

F [B^^) _ e"'(^o^^^")/^ R{B^^) _ e-^C^o^.^o"') 

The terms at the right side of the two equations in (4) 
are colliding systems independence and the two equa- 
tions are satisfied for any fixed B^ and 5^, respectively, 
which indicate that the variables of the two observables 
F {Bq) and R{Bq) can be separated as colliding system 
dependent term multiplied by scale impact parameter de- 
pendent term as: 

R{A,Eiab,Bl^) = ^''(A^iafc) X77«(6^) (5) 

F {A, Flab, B^) = e{A, Flab) X 77^(6^) (6) 

With the variable separable property, which is nontrivial 
and not common for all the experimental observables. 
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the excitation properties of the correlative observables 
for any scale impact parameter are the same. For better 
understanding the fitting parameters L and m, as well as 
the normalization factor pmbi the equ. (5) and equ. (6) 
are introduced into equ. (2), then we obtain: 



ln[ 



Pmb 



L X ln[ 



Pmb 

-Lxln[r;^(6o^)]-ln[,7^(6«)] (7) 



As for various colliding systems {A, Eiat) and scale im- 
pact parameters {b^,b^), the equ. (7) are always satis- 
fied, we have: 



In 



i'^iA.Elab 



L X ln[ 



e{A,Eiab), 



(8) 



Pmb Pmb 

m{b^, 6o«) = ln[r;^(6o^)] - L x \n[rj^{b^)] (9) 

The equ. (8) demonstrates that the colliding system de- 
pendent terms of the two correlative observables have 
been connected by introducing a proper normalization 
factor Pmb, which is also colliding system dependent and 
may be not unique for present correlation or even not 
necessary for other correlative analysis. In equ. (8), the 
universal fitting parameter L is uniquely determined by 
the colliding system dependent properties of the two ob- 
servables, that is the reason the L is not related to the 
bg and b^ for a given nuclear EoS (See Fig. 4). Thus, 
it is supposed to be a constant characteristic parame- 
ter to characterize the nuclear intrinsic properties. The 
variable separable properties of the two correlative ob- 
servables in equ. (5) and (6), and the analytic relation 
between colliding system dependent terms in equ. (8) are 
the origin of the correlation presented by equ. (2). As 
a consequence of the correlation for various scale impact 
parameters, the intercept parameter m can be expressed 
as equ. (9), which is the combination of the scale im- 
pact parameter dependent terms of the two correlative 
observables and without the cross terms. Derived from 
equ. (5), (6) and (9), the differential of the parameter 
m{bQ , 6^) can be written as: 



dm 
db§ 

1 dm 
Ldbl 



91n[r;«(6^)] _ dHR{A, Eiab,b^)] 



db^ 



db^ 



d\n[r,P{b^)] d\n[F{A,E,ab,bl^)] 



db^ 



db^ 



(10) 



(11) 



The differential of m{b^ ,b^) in equ. (10) and (11) are 
only related to b^ and 6^, respectively (See Fig. 4) and 
uniquely determined by the corresponding experimental 
observable. The colliding system dependent properties of 
the differential of experimental observables in equ. (10) 
and (11) can be used to validate whether the observables 
are variable separable or not, which is a necessary and 



not sufficient condition for the present correlation. Be- 
cause the equ. (2) is only the fitting equation for the 
two correlative observables and not all of the colliding 
systems exactly satisfy the equ. (2), the differential of m 
may be weakly dependent on specific colliding system. 

We have performed correlative analysis between the 
nuclear stopping and scale invariant nucleon sideward 
flow for various colliding systems and scale impact pa- 
rameters. The complication of colliding system depen- 
dence of single observable can be separated from the scale 
impact parameter dependence, and has been reduced to 
two universal fitting parameters L and to, which can be 
used to determine the nuclear EoS or other intrinsic prop- 
erties. The essential of the universal correlation behav- 
ior between nuclear stopping and scale invariant nucleon 
sideward flow from SIS to SPS energies may result from 
the pressure of the matter in HICs, which is dominated 
by nuclear EoS, in- medium NN cross section etc. and 
intimately connected to the nuclear stopping and nucleon 
sideward flow0, [H, [H, [13, [SI, [H . The strong corre- 
lation may indicate that the pressure produced in HICs 
may be also with the variable separable property as in 
equ. (8) and (9), and it may be broken down by the 
sudden change of the nuclear bulk properties, such as 
phase transition etc. As the phase transition mechanism 
is not explicitly involved in UrQMD model and it is also 
found in Ref. [1^ that the collapse of excitation func- 
tion of the sideward flow and elliptic flow can be used 
to probe the first order QGP phase transition. Thus, 
for qualitative analysis, it is predicted that the universal 
correlation may be broken down at different sets of beam 
energies for various colliding nuclei, which could serve as 
signals for the first order QGP phase transition. Fur- 
thermore, if the universal correlation pattern restores at 
much higher energies, where the crossover from hadronic 
phase to partonic phase would happen, the location of 
the critical point can be unitedly restricted by the map- 
ping of different colliding nuclei with the corresponding 
lower limit of the restoration energies on the QCD phase 
diagram. The real experimental data as well as the theo- 
retical calculation are expected to be compared with the 
UrQMD simulation results and the detailed correlation 
mechanism of the two experimental observables should 
be further studied. 
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